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We consider a double quantum dot coupled to two normal leads and one superconducting lead, 
modeling the Cooper pair beam splitter studied in two recent experiments. Starting from a micro- 
scopic Hamiltonian we derive a general expression for the branching current and the noise crossed 
correlations in terms of single and two-particle Green's function of the dot electrons. We then study 
numerically how these quantities depend on the energy configuration of the dots and the presence 
of direct tunneling between them, isolating the various processes which come into play. In absence 
of direct tunneling, the antisymmetric case (the two levels have opposite energies with respect to 
the superconducting chemical potential) optimizes the Crossed Andreev Reflection (CAR) process 
while the symmetric case (the two levels have the same energies) favors the Elastic Cotunneling 
(EC) process. Switching on the direct tunneling tends to suppress the CAR process, leading to 
negative noise crossed correlations over the whole voltage range for large enough direct tunneling. 

PACS numbers: 73.23.-b, 74.45.+C, 72.50.+m 

I. INTRODUCTION 

At low temperatures, electron flow in mesoscopic systems bears analogies with the propagation of photons in 
quantum optics devices. The fermionic analog of the Hanbury Brown and Twiss experiment for photon is an example: 
negative current-current correlations demonstrate that the statistics of current carriers in microstuctures corresponds 
to a degenerate Fermi gas. Over the last two decades, the issue whether all fermionic systems should demonstrate 
antibunching has been addressed. It was predicted that if electrons are injected from a superconducting lead, into a 
fork consisting of two normal metal leads, then positive current-current crossed correlations could be observed^—. This 
phenomenon, called Crossed Andreev Reflection (CAR), has since been interpreted as originating from the splitting of 
the constituent electrons of a Cooper pair from the superconductor into the two normal leads. Because this Cooper pair 
is a singlet pair of electrons, the two electrons which are injected into these leads should preserve their singlet nature, 
therefore providing a modern example of the non-local character of quantum mechanics in nanoelectronics£~— . Bell 
inequality measurements based on noise crossed correlations have been proposed to detect this entanglement £iLr— 




FIG. 1: (color online) Double quantum dot coupled to normal/superconducting leads. 



Over the last decade, CAR has been studied extensively for a variety of geometries. The early proposals pointed 
out that positive crossed correlations can be generated either by selecting the energies of the outgoing electrons, or 
alternatively by filtering their spins.— Some theoretical works consider a superconductor connected in two separate 
locations to normal metal leadsJ^— Such works address the effect of the separation between the normal metal 
leads on the CAR signal, and they specify which voltage configuration results in positive crossed correlations. Other 
theoretical works consider hybrid structures where the normal leads are taken to be ferromagnets or half metals&22r— , 
semiconductors 25 and Luttinger liquids ! 26 ' 27 Recently the effect of Coulomb interaction has been studied in some of 
these systems ! 28 ' 29 
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On the experimental side, noise experiments with normal superconducting devices constitute a real challenge. A 
successfull experiment demonstrated that the Fano factor of a single NS junction is 2, corresponding to the charge of 
a Cooper pair.— One of the intrinsic difficulties lies in achieving controlled, good quality contacts between the two 
arms and the superconductor It has also been suggested that non-local effects in the current (Andreev drag) could 
be probed by placing two separate contacts to the superconductor. The first drag experiments were performed in this 
geometr y 32 ! 33 , but one true challenge remains that the signal for drag effects is very weak, and again that symmetric 
contacts are difficult to achieved 

Two recent experiment; ; 35 ! 36 have provided evidence of Cooper pair splitting in devices consisting of a supercon- 
ducting finger placed on the bulk of a single nanowire, whose both ends are connected to metallic leads (see Fig. f). 
Because both the superconducting and the metallic leads are placed on top of the nanowire, two quantum dots are 
generated on both sides of the superconductor, and their energy levels can be controlled with the help of additional 
gates. This system constitutes a tunable Cooper pair beam splitter with a relatively good degree of symmetry. Dif- 
ferential conductance measurements showed appreciable non local signal which could be attributed to CAR. Some 
results also suggest that local interaction and direct tunneling between the dots may play an important role. In this 
work, we choose specifically to focus on the geometry of Refs. [Ill and The main goal of the present paper is to 
compute the branching currents of this device and their crossed correlations. These quantities depend strongly on 
the energy configuration of the dots as well as the presence of direct tunneling between them. In particular we study 
which parameters have to be optimized in order to achieve Cooper pair splitting. 

While exploring the above physics, we put some emphasis in presenting the formal aspects of the calculation, which 
are achieved starting from a microscopic Hamiltonian describing a BCS superconductor and normal leads, all coupled 
via tight binding tunnel Hamiltonian to the dots. Using a path integral formalism, the leads degrees of freedom 
are integrated out, which is equivalent to resumming the Dyson series in a perturbative Green's function approach. 
This allows to obtain a rather general formula for the current in terms of the Green's function of the dots. More 
importantly, we are also able to derive a general expression for the noise crossed correlations, in terms of a two particle 
dot Green's function. 

The outline of the paper is as follows. In Sec. UH we introduce the model for the Cooper pair splitter. The 
currents and the crossed correlations are computed in Sec. IIII1 and finally we apply our results to various situations 
(symmetric/antisymmetric case and with/without direct tunneling) in Sec. IIVI to describe in which regime Elastic 
Cotunneling (EC) or CAR are optimized. We conclude in Sec. [V] 



In this section we introduce the Hamiltonian of the double dot system and derive the expression for the tunneling 
self-energy. To do so, and in order to clarify notations, we start by studying the single dot formalism^ - — 



We consider a quantum dot with a single level e, which is coupled via tunneling amplitudes tj (J — L,R,S) to 
two normal leads and one superconducting lead with superconducting gap A (see Fig. We label the applied 
bias voltage on the left (right) lead Vt (Vr), where Vj is measured with respect to the chemical potential of the 
superconducting lead whose voltage is set at Vs = 0. For simplicity we work with physical dimensions corresponding 
to h = I and e = I. The Hamiltonian of the total system reads 



II. MODEL 



A. Single dot formalism 



1. 



Hamiltonian 




(1) 



.7 



where the dot Hamiltonian is given by 




(2) 



<7=t,4- 



and the lead Hamiltonians are expressed in terms of Nambu spinors 




(3) 



k 
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FIG. 2: Single quantum dot coupled to normal/superconducting leads where t L / R (ts) is respectively the tunneling amplitude 
between the dot and the normal (superconducting) lead. 



with a z , a x Pauli matrices in Nambu space and 

In the case of normal leads, Aj is zero. The tunneling Hamiltonian, which is responsible for transfer of electrons 
between the leads and the dot reads 

H T (t) = J^*] k T j (t)d + h.c, (5) 

jk 

where we introduce the Nambu spinor of the dot electrons 
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We include the voltage dependence in the tunneling term using the Peierls substitution (performing a Gauge trans- 
formation in order to represent the bias potential as a vector potential) . The tunneling amplitude thus reads 

T j (t)=t J a z e Ur 'f v > dt . (7) 
We now introduce the bare Green's functions of the dot (in the absence of tunneling) 

G s s '(t,t') = -i(T c {d s (t)d<t s '(t')}) Q , (8) 

where Tq is the time ordering operator along the Keldysh contour, s, s' labels the position of the times on this 
contour. The quantum mechanical averaging is performed with respect to the Hamiltonian without tunneling (. . . ) = 
ZQ lr Tr{e- 0H ° ...} with Z = Tr{e~^ H °} and H = H D + £\ H r The Green's function dressed by tunneling is 
written as 

G ss '(t,t>) = -i (t c {s(ooK(t)rft'y)}) , (9) 
where S(oo) is the evolution operator along the contour 



\-i [ + dt ]T 

I J -°° s=+,- 



S(oo) = T c cxp{ -i I dl ) rril^(l) } . (10) 
and t z is the z Pauli matrix in the Keldysh space. 
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2. Averaging over the leads 

In this section we calculate the self-energy associated with the tunneling between the dot and the leads. Because the 
lead degrees of freedom appear quadratically in the total Hamiltonian, the evolution operator can straightforwardly 
be averaged over such leads 



( S (°°)) leads =T C exp 

where we introduce the spinor in the Nambu-Keldysh space 



i /<fti<ft 2 #(ti)iMti,ta)<i(t2) 
Jc 



*» - ( tt ) ■ 

Here the self-energy associated with the tunneling takes the form 

£ T (ti,t 2 )= J2 %{tiM)= tf{t x )®T z )g j {ti-t 2 )(T z ®T j {t2)) , 

j=L,R,S j=L,R,S 

where Sy, gj are matrices in Nambu-Keldysh space, and gjit — t') — —i J2k i^c {^jk(*)^jfc(*0 }} * s tne l° ca l Gr 
function of electrons on lead j. In the litterature^, they are typically given in the rotated Keldysh space, 

1 f 1 



(ii) 



(12) 



(13) 



%RAK _ 



Lf z gL 1 , with L = 



v^V 1 1 



(14) 



where f z = t z ® 1, and 1 is the unity matrix in the Nambu space. With this rotation the Green's function matrix 
can be simply expressed in terms of the advanced, retarded and Keldysh components 
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RAK 



9 R 9 K 



g A 

The components of the leads Green's functions in the rotated base have the following form 



(15) 



g. (w) = 7n/(0) u-a ■ ■ 



, uj-l-\-Aj-a x 



5 fM=(l-2/ u )( fl f(,)- 5 /( W )) 



(16) 



± sign(w) A u) 



where C 



R,A 



<JA) - a; 2 , 



M > A, 



M < Aj 



(17) 



with v(fi) the density of states of the normal lead at the Fermi level and fa, the Fermi distribution. We rotate back 
in Keldysh space the matrix for the lead Green's function and we get the following formula for the self-energy £j in 
time-domain 



E J -(ti s t 2 ) = r j 



duj 
2^ 



e -iu.(ti-t 2 ) e -«r.VGifi[ w . 1 - Aj ■ a x ]e +ia ' v ^ 2 



(18) 



e(A,-|^h v©(|w| - Aj) / 2/ w - 1 -2/, 



A|-c 2 



-r z + zsign(w)- 



- A 2 



-2/- w 2/ w - 1 



where Tj = ttv(0) is the tunneling rate between the dot and the lead j. 
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FIG. 3: Double quantum dot coupled to normal/superconducting leads. 



B. Double dot formalism 



For the remainder of this paper we focus on a system of two single-level quantum dots coupled to one superconducting 
lead and two normal leads (see Fig. [3]). Such a system was studied experimentally in Refs. |H and HI, using a carbon 
nanotube or a nanowire attached at both ends to normal metal leads with a superconducting electrode in the middle. 

The double dot formalism is an extension of the single dot one where there appears a new matrix structure in Dot 
space. The Hamiltonian of the total system reads 



H = H Dl + H D2 + H DlD2 +J2 H i + H Ti (*) + h t 2 (*), (19) 

i 

where Hn x and Hr> 2 are the Hamiltonians of the two dots 

H Da = e a d Ld**, (20) 

<T=f4 

and the tunneling between the dots is conveniently written using Nambu spinors and the tunneling amplitude td 

Hd 1 d 2 = tdd\(J z d2 + h.c. . (21) 
These three terms can easily be combined under the following form 

H D = <? ^ ® a J, (22) 
where t? = (d\ d^j is a spinor in the Nambu-Dot space. The tunneling between the dot a and the leads reads 

Hr a (t) = T ^ W da + h - c - ' ( 23 ) 

jk 

where Tj a {t) — tj a a z e lcr * / ^ dt ; ancl tj a corresponds to the tunneling amplitude between the dot a and the lead j. 
This Hamiltonian corresponds to a system where every dot is coupled to every lead. In order to reproduce the system 
of Fig. [3J it is necessary to set to zero the tunneling amplitudes t^2 and tux coupling the dot 2 with the left lead and 
the dot 1 with the right lead. 

Now the Green's function of the dot electrons has a new form 

G s a Ut,t') = ~i (Tc {d£(t)d#(i')}) , (24) 
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with a, a' corresponding to the dot index. As in the previous section the self-energy associated with the tunneling 
between the dots and the leads is calculated by averaging the evolution operator over the leads degrees of freedom 



( S (°°)) leads =TceX P 



i / dhdt 2 S(h)£ T (t u t 2 )d{t 2 ) 



— oo 



(25) 



d-(t) 

t (t -t) = V ( C 7 ?!^)® 7 "*)^'^ -ta)fo®75i(i a )) {T^{t 1 )®T z )g ] {t 1 -t 2 ){T z ®r ]2 {t2))\ . . 



where d(t) = I I is a spinor in the Nambu-Dot-Keldysh space and the self-energy is given by 



j 



where g.j and St are matrices in Nambu-Keldysh and Nambu-Dot-Keldysh space respectively. Each element £ Jce ^(ii — 
t 2 ) of the self-energy matrix (|26p in the Dot space can be obtained from Eq. (fT5|) by replacing Tj with Tj a p = 

III. CURRENTS AND CROSSED CORRELATIONS 

In this section we derive the currents between the two dots and the various leads as well as their crossed correlations 
using the dot electrons Green's function and the tunneling self-energy calculated previously. 

A. Currents 

The current from the dot a into the lead j reads 

W*) =i$^¥} fc ff,7} a (f)d Q + h.c. . (27) 

k 

The average current does not depend on which branch of the Keldysh contour the time is chosen, and therefore can 
be expressed as (Ij a ) = (Ija{t + ) + Ij a (t~)) /2 where is the time on the upper/lower branch of the contour. In 
order to compute these it is convenient to introduce counting fields T]j a {t) which appear in the tunneling amplitudes 
as Tj a (t) — ¥ Tja(t)e , ' T *®' 7zr,ja< ' t ^ 2 ■ The average current from the dot a into the lead j can then be calculated as the 
first derivative of the Keldysh partition function 



Z[0] 6 Vja (t) 



(28) 

?7=o 



where Z[rj\ = (S(oo,r])) and 5(00,7/) is the evolution operator in which the counting fields were introduced. After 
performing the derivative we obtain the following result for the current 

(I ja ) =^Tr|(r z ® a z )£°°dt> (G(t, t')t 3 {t' , t) - £j(t,t')G(t> \t)) \ , (29) 

where "Tr" corresponds to the trace in Nambu-Keldysh space. Going to the rotated Keldysh space, the average 
current can be re-expressed in terms of the advanced, retarded and Keldysh components as 

(Ija) =^tr y + °° dt' (G R (t, t')tf(t', t) + G K (t, t')tf{t', t) - tf(t, t')G K (t', t) - Ef (t, t')G A (t', t)) ^| , (30) 

where "tr" corresponds to the trace in Nambu space. While the self-energy Sj can be obtained from the results of 
the previous section, the Green's function G remains to be determined. To do this we write the Dyson equation in 
the frequency domain and we obtain the various components of G in the rotated Keldysh space 

G^iuy 1 = G^iuj)- 1 - S^H, (31) 
G K (u) = Gff(u) + G R {u)t% {uj)G a {uj), (32) 
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with 



r,R/A, v_i _ fut - e%a z -t d cTz 
G* M - 0. 

In terms of these Fourier transformed functions, the current can be rewritten as 

( +00 

(lia) = tr I a z [ ^Re [(G R (u)tf (u) + G K (uj)tf(u J ) 



(33) 
(34) 

(35) 



where we used that the Keldysh components are anti-Hcrmitian while the advanced and retarded components are the 
Hermitian conjugate of one another. We can now calculate the current using the Dyson equations (|3T |) - (f32|) and the 
components of the tunneling self-energy obtained in the previous section which read in the frequency domain 



= ±i 



. f 



r 



R22 



'1 , 



= ±i 



.■fr Ln 







(w) = X* /R (u) 



Ts21 T522 



1 



Ef(w) = -2i 





r^22 



tanh (^_v«)) 







t a nhf^±^ 



1 -a 



rsil T5i2 
r521 Tg 2 2 



(36) 
(37) 
(38) 

(39) 

(40) 
(41) 



where we focused on the case Vs — allowing us to simplify the expression for the self-energy Eg, and introduced 



x A/R (u) = Q ( A -M)" ± , e(^|-A) H 



VA 2 - w 2 



Vw 2 - A 1 



(42) 



Xf (w) = -2i 



e(M-A) h 
V^ 2 - a 2 V 2 



(43) 



B. Crossed correlations 



We follow a similar approach to the one developed to calculate the currents. Here however the crossed correlations 
involve two current operators evaluated at different times which in general do not commute. Usually the measurement 
procedure dictates which combination of the current-current correlators is involved in the expression of the measured 
noise.— For this reason, we need to compute the unsymmetrized correlator, which can subsequently be manipulated 
to be applied to a desired measurement procedure. 

We introduce a partition function Z[rj\ = (S(oo,r))) depending on a counting field rjj as (t) where j = L,R,S 
corresponds to the different leads, a = 1, 2 to the two dots and s = ± to the branches of the Keldysh contour. This 
new counting field enters in the tunneling amplitudes as 

T ja ->T ja eE. i **®°'*i<»W, (44) 

where we introduced 7r-matrices in Keldysh space 



n+ = (0 o) ,7r - = (0 -l) 



(45) 
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Now the crossed correlation can be calculated as a second derivative of this partition function over the counting fields 

f 6 2 Z[rj\ 

Performing the average over the leads allows us to write the current-current correlation function 



/ ~ Z[0] 5r, ia _(t)Sr, J0+ (t>) 



(46) 



1^)1^')) = J dt ^ E E ^ 



X ^ 1 jMl)Si; J ,(t',t 2 )^ + (*l,t a ,t,*') -^ < T^ oa (*,*l)£^ ^ V a ai( t2 ' t, )^?' 



(Ti (T 2 0"2 0"i O"! £72 CT 2 / 

where 



te s 4%(ii,i2,t3,t4) - - (T C {< CTl (^)C CT2 (i 2 )fefe)4 s 4 4 . 4 (i4)}) (48) 

CTiO-2CT3(T4 

is the two particle Green's function of the dots electrons and Ej 1 *^^^ (ti,t 2 ) is the matrix element of the tunneling 
self-energy associated with the lead j. In the general case where Coulomb interaction is present on the dots, the 
two-particle Green's function can be expressed in terms of the dressed single particle Green's function and the full 
vertex function. However in our non-interacting case the two-particle Green's function reduces to 

Ka 1 a 2 %a 4 (h,t 2 ,t 3 ,t i ) = G^a^ai (*i , U)G a 2 %v 2<T3 (*2 , *3 ) ~ G^ % ai cr 3 (* 1 ^G^^ <Jl (t 2 , U) , (49) 

fJlfT2^3^4 

where G* 1 i s J 2Cricr2 {t\, t 2 ) is the a±a 2 matrix element of the single particle Green's function introduced in Eq. (|24l) . 
Substituting Eq. (1491) into Eq. (1471) we obtain the irreducible part of the current-current correlator 

SiajfS&t') = (I ia (t)I jp (t')) - (I ia (t)) (I^(t')) 

+ 00 

dhdhTr { (tt_ <g> o^X^t, ti)6(ti, t'))a/9(7r+ ® 0-«)(Ej (*', t 2 )G(t 2 ,t)) 0a 

-OO 

+ (tt_ <g> a z )(G(t, t 1 )t j (t 1 ,t')) al3 (ic + ® ^)(G(t', i 2 )Si(« 2 , 

- (tt_ ® o-,)(Ei(i, ti)(5(ii, i 2 )Sj (i 2 , 0)o/9(t+ ® a z )G 0a (t', t) 

- {K-®a z )G aP {t,t')(K + ®a z ){t J {t'M)G{t 1 M)%(t2,t))p a } . (50) 

Performing a rotation in Keldysh space, the irreducible part of the current-current correlation function can be ex- 
pressed in terms of the advanced, retarded and Keldysh components as 

+00 +00 

Si a ,jfi(t,t') = -~Re J J dt x dt 2 

— OO — OO 

x tr L z ( Ef G A + Ef G K - Ef G A + t?G a ) ^ a z ( Ef G A + Ef G A ' + Ef G A - Sf G R Y" 

IV* / (t,ti)o(ti,t') V J J J J / (t',t 2 )o(t 2 ,t) 

-a z ( Ef G*Ef + Sf G A Ef + Ef G K Ef - Ef G A t A + ± R G R ±f) ^ a z (c K +G A ~ G R ) "° \ . 

Finally we take the Fourier transform of this expression and we obtain the irreducible part of the frequency dependent 
current-current correlation 

+00 

SiaM") = ~l Re J ^ {^?G A + ^fG K ~ ^ A + SfG R )^ a z (tjG A + ZfG K + t A G A - EfG*)^ 

— OO 

-a z (EfG fi Ef + Ef G A tf + t R G K tf - Ef G A tf + Ef G R tf )^ a* (g x + G A - G R )^ | . (51) 

In the following section we compute numerically the current (|35[) and the crossed correlation (|51[) for various regimes 
and we comment on the results. 
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IV. RESULTS AND DISCUSSION 
A. Transport without tunneling between the dots 

To begin we present the three dominant processes which can occur in such a system. Then we study two situations: 
the anti-symmetric case when the energy of the two dots have the opposite position (with respect to the chemical 
potential of the superconducting lead at Vs = 0) and the symmetric case when the energy levels of the two dots are 
the same. All energy scales in this section are in units of A and in the subgap regime. In all the following results we 
focus on the low temperature regime (/? ^> 1/A) as the only effect of temperature is to smooth the signal. 

1. Dominant electron transfer processes 




(a) Direct Andreev Reflection {b) Crossed Andreev Reflection (c) Elastic Cotunneling 



FIG. 4: Dominant electron transfer processes 

Direct Andreev Reflection (DAR) (see Fig. [4^,) This process involves an electron which is incident from a 
normal lead at energies less than the superconducting energy gap. The incident electron forms a Cooper pair in the 
superconductor with the retroreflection of a hole of opposite spin and momentum to this incident electron. This 
is equivalent to two electrons with opposite spin and momentum from the same lead forming a Cooper pair in the 
superconductor. The inverse process corresponds to the destruction of a Cooper pair in the superconductor with two 
electrons propagating in the same lead. 

Crossed Andreev Reflection (CAR) (see Fig. [4J5) Crossed Andreev reflection occurs when two spatially 
separated normal leads form two separate junctions with a superconductor (the separation between the two normal 
leads is assumed to be smaller than the BCS superconducting coherence length). In such a device, retroreflection 
of the hole from an Andreev reflection process, resulting from an incident electron at energies less than the 
superconducting gap at one lead, occurs in the second spatially separated normal lead with the same charge transfer 
as in a normal Andreev Reflection (AR) process to a Cooper pair in the superconductor. For CAR to occur, electrons 
of opposite spin and energies must exist at each normal lead (so as to form the pair in the superconductor). The 
inverse process corresponds to the destruction of a Cooper pair with two electrons propagating in opposite leads. 

Elastic Cotunneling (EC) (see Fig. Hfc) Elastic cotunelling is the quantum mechanical tunneling of electrons 
between the normal leads via an intermediate state in the superconductor. It does not lead to Cooper pair creation 
or annihilation in the superconductor. 

According to the configuration of our system one or several of these processes occur. The goal is to study which 
configuration facilitates which process. 

2. Anti-symmetric case 

In this section we focus on the geometry where the two dots levels are opposite with respect to the superconducting 
chemical potential. We fix the voltage of the right lead and we vary the voltage of the left one (Fig. [5]). 

In Fig. [5] we plot the current which flows in the two normal leads (Ili is the current between the dot 1 and the 
left normal lead and Ir2 is the current between the dot 2 and the right normal lead) . We see that for voltages below 
ei (the energy level of the first dot) both currents have the same sign (with our notation currents which enter the 
leads are positive) and essentially the same amplitude, (strictly speaking, we observe a small deviation between the 
two currents around Vl = £2)- Above t\ the currents have a reduced but comparable amplitude but their signs are 
opposite. 
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FIG. 5: Anti-symmmetric case: the energy levels of the dots are opposite with respect to the chemical potential of the 
superconducting lead whose voltage is set at Vs = 0. The voltage of the right lead Vr is fixed and we vary the voltage of the 
left one Vl- 
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FIG. 6: (Color online) Currents (arbitrary units) as a function of the voltage of the left normal lead Vl for ei = 0.5, ei = —0.5, 
Vr = -0.7, /3 = 100, t L i = t R2 = tsi = t S 2 = 0.2 and t L 2 = tm = 0. 



When the voltage of the left lead is smaller than 62, the dominant process is the Crossed Andreev Reflection (CAR) 
in the sense that a Cooper pair from the superconducting lead is split and its constituent electrons are injected in the 
two leads. Strictly speaking, below e 2 Direct Andreev Reflection (DAR) gives also a contribution, albeit a minor one. 
This can be understood by looking at the density of states of the dots: 

p a = -Im(G A ) ttctct , (52) 

where a = 1, 2. p\ is plotted in Fig. [7] It contains a sharp double peak at ei, and a much smaller peak at —e\ = e 2 
which originates from the proximity effect of the superconducting lead. In this situation p 2 (not shown) is exactly the 
symmetric of p\ with respect to the uj — axis. As electrons can tunnel through these two resonances in the same 
lead, this explains the presence of a weak DAR process. 

The CAR and DAR processes thus give currents (Ili and Irz) with the same sign because electrons injected from 
the superconductor only end up in the two normal leads. Elastic Cotunneling (EC) does not give any contribution 
because the voltages of each normal leads are not large enough to allow quasiparticle injection in the superconductor. 
Note that the amplitude of the generated currents by CAR is dominant because this configuration is optimal for the 
process of resonant electron transfer through the (large) peaks of the density of states p\ and p 2 (at opposite energies). 

As soon as V L > — ei, the two currents deviate slightly from each other because EC comes into play. Indeed, an 
electron from the left lead can tunnel into the superconductor via the (small) resonance in p\. It may continue its 
way to dot 2 using the (large) resonance in pi which is located at the same energy. Note at the same time that DAR 
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processes continue to operate, but only for the right lead. This explains why Ili > Ir2 in the region [— ei, ej.]. 

When the voltage of the left lead is larger than ei , the CAR process is strongly suppressed because Vl is placed above 
all resonant levels of p\ and p2- Electrons injected from the superconductor are prohibited from entering the normal 
left lead. The only allowed processes are thus EC and DAR. DAR process from the left lead injects Cooper pairs in 
the superconductor, while they destroy Cooper pairs which are injected as electrons in the right lead. In addition, EC 
contributes to transfer electrons from the left lead to the right lead via the superconductor. The amplitude of such 
currents is reduced (compared to the CAR amplitude in the interval [— e\, t\\) because both DAR and EC require one 
electron transfer through a "small" resonance of p\ and/or pi- 
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FIG. 8: Current-current crossed correlation Sli,r2(0) (between Ili and Ir2) at zero frequency (arbitrary units) as a function 
of the voltage of the left lead Vl for the same setup as in Fig. [6] 

In order to confirm these observations, we plot in Fig. [8] the current-current crossed correlation at zero frequency 
between Ili and Ib.2- The crossed correlation is positive until t\ where the signal drops to zero and becomes eventually 
negative. Strictly speaking there is a small structure at — a due to the (small) resonance of the density of states. At 
ei the cross correlation has a small secondary peak which originates from the (large) double peak resonance in the 
dot density of states. For voltages larger than ei, the negative signal has a much weaker amplitude than the positive 
one (below ei) because the associated currents are smaller. 

We now give a physical interpretation of these results. When the voltage is smaller than ei, the dominant process 
is CAR and the crossed correlations are positive because the two electrons of the same Cooper pair are split and end 
up in opposite leads. When the voltage is larger than ei the dominant process is EC and the crossed correlations are 
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negative. An electron injected from the left lead, which tunnels trough the superconductor then enters the right lead. 
Note that the generated currents by DAR do not contribute to the crossed correlation because the two Cooper pairs 
which are injected in each side of the superconductor are independent. 



3. Symmetric case 

We now focus on the second geometry where the two dots levels have the same position (see Fig. HJ). Again we fix 
the voltage of the right lead Vr and we vary the voltage of the left one Vl- 




FIG. 9: Symmetric case: energy level of the dots are the same. The voltage of the right lead Vr is fixed and we vary the voltage 
of the left one Vl- 
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FIG. 10: (Color online) Currents (arbitrary units) as a function of the voltage of the left normal lead Vl for e\ — 0.5, ei = 0.5, 
V R = -0.7, /3 = 100, t L i = t R2 = tsi = t S 2 = 0.2 and t L2 = tm = 0. 

In Fig. [10] we plot the currents in the two normal leads. We see that contrary to the previous case, the two currents 
have a reduced amplitude below Vl — ei and this amplitude increases above this voltage. To be more precise, the 
two currents are negative and identical below Vl = — £i- In the interval [— ei,ei] the currents remain negative but 
they differ slightly. Above Vl = ei, these currents bear opposite sign and their amplitude has increased. 

When the voltage of the left lead is smaller than — ei, the dominant processes are CAR and DAR. These two 
processes distribute electrons in the two leads in an equivalent manner, hence they are equal. The smallness of their 
amplitude can be explained as above from density of states considerations (see Fig. [TTj) . Note that in this geometry 
Pi = P2- both contain a large peak at t\ and a much smaller one at —t\. For both CAR and DAR processes, an 
electron transfer through the small resonance is required, which explains their amplitude. In the interval [— ej.,ei], 
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DAR processes from the left lead are suppressed altogether, but they still operate with the right lead. The CAR 
process also contributes (the two currents still bear the same sign) , and it involves electron transfer through both the 
"large" resonance of dot 1 and "small" resonance of dot 2 (hence the reduced amplitude) . Note that in addition there 
arises a small contribution from EC processes, but this contribution is less important because the electron has to be 
transfered through both "small" resonances of the dots. 

When the voltage of the left lead is larger than t\ , the dominant process is EC because electrons of the left lead can 
be injected in the first dot (through its "large" resonance) they tunnel trough the superconductor and then they end 
up in the right lead via the second dot (through its "large" resonance also) . In this configuration the DAR is always 
present but its amplitude is reduced, once again because of density of states arguments (passage through a "small" 
resonance). These two processes give currents with opposite sign because either electrons from the left lead go into 
the right lead (EC processes), or two electrons from the left lead form a Cooper pair in the superconductor (DAR 
from left lead), or Cooper pairs are injected from the superconductor into the right lead (DAR into right lead). CAR 
is killed because injection of electrons in the left lead is prohibited since the voltage is larger than the dot 1 energy. 
The amplitude of the generated currents in this configuration is important because this corresponds to the optimal 
configuration for EC. Electron transfer then exploits the maximum of the density of states for both dots. 
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FIG. 12: Current-current crossed correlation Sli,r2(0) (between Ili and Ir2) at zero frequency (arbitrary units) as a function 
of the voltage of the left lead Vl for the same setup as in Fig. 1101 

In Fig. [12] we plot the current-current crossed correlation at zero frequency between Ili and Ir2- The correlations 
are positive below ei then become negative, with a large amplitude. Below ei the amplitude of the signal is reduced 
and it has a structure at Vl — — ei- 
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When the voltage is smaller than t\ , the dominant process is CAR thus the crossed correlations are positive and the 
explanation is the same as in the anti-symmetric case. The reduced amplitude of the crossed correlation is explained 
by the fact that one electron has to be transfered through a "small" resonance. The feature at Vl = — ei corresponds 
to the onset for EC process through the small resonances of the two dots. Above Vl — t\ electrons can pass through 
both large resonances of the dot which explains the large (negative) amplitude of the signal. The generated currents 
by DAR don't contribute to the crossed correlation because the two Cooper pairs which are injected in each side of 
the superconductor are independent. 

To summarize we can facilitate the EC regime if we have the same energy level position for the two dots (Symmetric 
case) and the CAR regime is facilitated if we have opposite energy level for the two dots (Anti-symmetric case). In 
the following section we are going to be interested in the effect of direct tunneling between the two dots, which is 
relevant in recent experiments , 35 i 36 

B. Transport with tunneling between the dots 

In this section we allow tunneling between the two dots (t<j ^ 0) and we study the effect of this manipulation on 
the currents and the crossed correlation at zero frequency. 

1. Antisymmetric case 

The main effect associated with tunneling between the dots is that it modifies the density of states of each dot. In 
Fig. [T3]we see that for the anti-symmetric case, the weight of the small resonance (at — e x for dot 1, and at ei for 
dot 2) is increased by such tunneling. The large resonance (at ei for dot 1, and at —e% for dot 2) now acquires an 
asymmetric double peak structure, for positive energies, which disappears for sufficiently strong tunneling between 
the dots. 
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FIG. 13: (color online) Density of states of the dot 1 (arbitrary units) for various values of the tunneling between the two dots 
and for ei = 0.5, £2 = —0.5, thi = tm. = tsi = ts2 = 0.2 and thi = tm = 0. 

The two currents are plotted in Fig. [14] For Vl < — ei, both currents are negative with a large amplitude associated 
to CAR, as in section HV A 21 Contrary to section HV Al they have a different amplitude (Ir2 > Ili)- There are two 
dominant processes into play. On the one hand, the presence of tunneling between the dots allows to transfer electrons 
from right to left lead without passing through the superconductor. On the other hand, as the (small) resonance of 
dot 1 has been increased, there is also now the possibility for EC from dot 2 to dot 1 (this effect was not noticeable 
in the absence of tunneling between the dots). In addition, there is also the possibility for DAR process which injects 
electrons in both leads. 

Increasing Vl beyond — ei, CAR is still dominant, but the currents cross (Ili > Ir2 as in section BV A 2|) because 
both direct tunneling and EC from dot 1 to dot 2 come into play. The difference Ili — Ir2 is larger than in the 
absence of tunneling between dots because of the increased weight of the small resonance in the density of states. In 
this situation, there is also a contribution for DAR processes in the right lead. 




15 



Currents 



0.02 


- ** 
* 
* 
* 
* 
• 
> 


1.0 -0.5 


0.5 1.0 

■ 
■ 

:/ 

// 
y i 


/li 

Ir2 


-0.02 
-0.04 


* 

* 

"7."" 


......HOfi 




i 
i 
i 
i 







FIG. 14: (color online) Currents (arbitrary units) as a function of the voltage of the left normal lead Vl for ei = 0.5, e-2 = —0.5, 
Vr = —0.7, /3 — 100, tLi = tR2 = tsi = ts2 = 0.2, tL2 = tin — and the tunneling between the dots td = 0.2. 



Finally, for Vl > £i CAR is suppressed. In this regime, the main processes are direct tunneling and EC from dot 
1 to dot 2, as well as DAR from both the left lead to the superconductor, and the superconductor to the right lead. 
The difference Ili — Ir2 is further increased by both direct tunneling between dots and the increase of EC associated 
with density of states effects. 
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FIG. 15: (color online) Current-current crossed correlation Sli.r2 (0) (between Ili and 7_R2) at zero frequency (arbitrary units) 
as a function of the voltage of the left lead Vl for various values of the tunneling between the two dots and for ei = 0.5, 
<L2 = -0.5, Vr = -0.7, /3 = 100, t L i = t R 2 = tsi = t S2 = 0.2 and t L 2 = t m = 0. 

In Fig. 1151 we plot the zero frequency current-current crossed correlation as a function of the voltage of the left 
normal lead for various values of the tunneling between the two dots. For sufficiently weak tunneling parameters, the 
general tendency is to favor positive correlations for Vl < £i and negative cross correlations for Vl > £i, signaling a 
transition from CAR processes to EC processes. When the tunneling between dots is increased, the crossed correlations 
are shifted towards negative values, which results in a reduced positive signal in the CAR regime and an increased 
(negative) signal for the direct tunneling and EC regime. Beyond td — 0.7, the positive crossed correlations signal 
disappear completely. 

We further discuss the structure associated with the dot resonances at ±ei. Increasing the direct tunneling between 
the dot affects their density of states (see Fig. fl~3|) - As mentionned above, the "large" resonance at ±ei (for dot 
1/dot 2) has a double peak structure which disappears when td is increased, while the "small" resonance acquires a 
double peak structure when increasing td- This constitutes the justification for the side peak in the current-current 
correlation near Vl — ei to be smoothed out in the presence of strong tunneling. At the same time, for intermediate 
tunneling td = 0.2, the new double peak structure in the density of states of dot 1 at Vl = — £i creates a structure in 
the (positive) current crossed correlation. Further increasing td, a peak in the crossed correlation around Vl = — ci 
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is generated: this peak is shifted towards Vl < — ei in accordance with the shift in energy which is observed in the 
density of states of dot 1. This peak is the last stronghold for the observation of positive crossed correlations: for 
voltages below (above) it, direct tunneling between the dots transfers electrons from the right to the left lead (from 
the left to the right lead). We therefore interpret the presence of this peak as the point where the main process in 
competition with CAR, which is the direct tunneling, changes sign. 



2. Symmetric case 



The density of states (which is the same for dot 1 and dot 2) is drastically different in the case where the two dots 
have the same level position when the tunneling between dots is switched on (see Fig. I16|) . For moderate tunneling 
between the dots, it contains a (large) double peak at positive energies which is centered near e%, and a much smaller 
peak (not shown) close to —e\. Further increasing the tunneling, the double peak becomes well separated, with one 
side peak shifted towards negative energies, while the other side peak approaches the superconducting gap. 




FIG. 16: (color online) Density of states of the dot 1 (arbitrary units) for various values of the tunneling between the two dots 
and for ei = 0.5, £2 = 0.5, tz,i = tm = tsi = ts2 = 0.2 and tia = tm = 0. 
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FIG. 17: (color online) Currents (arbitrary units) as a function of the voltage of the left normal lead Vl for ei = 0.5, £2 = 0.5, 
Vr — —0.7, f3 = 100, tLi ~ tin = tsi ~ ts2 = 0.2, tL2 = im = and the tunneling between the dots td = 0.2. 

The two currents are plotted in Fig. [T7] For t& = 0.2 the currents In and Iri remain negative for Vl < 0, which is 
symptomatic of CAR and DAR processes, but their amplitude is drastically reduced because of the lack of weight in 
the density of states for negative energies. For Vl > 0, the two currents acquire an opposite sign and we observe two 
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steps which correspond to the double peak structure in the density of states. This corresponds to the situation where 
both direct tunneling and EC processes contribute in the same direction. Further increasing the coupling between 
the dots enhances these features. 
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FIG. 18: (color online) Current-current crossed correlation Sli,r2(0) (between Ili and Ir2) at zero frequency (arbitrary units) 
as a function of the voltage of the left lead Vl for various values of the tunneling between the two dots and for ei = 0.5, 
e 2 = 0.5, V R = -0.7, $ = 100, t L1 = t R2 = tsi = t S2 = 0.2 and t L2 = tm = 0. 



The crossed correlations are plotted in Fig. [18] for several values of the tunneling amplitude between the dots. The 
sign is positive but vanishingly small for Vl < 0, which is again the consequence of the drastically reduced density 
of states for negative energies. The cross correlations decrease monotonously when Vl is increased, but unlike the 
case of zero coupling between the dots the crossover from the CAR dominated regime to the EC dominated regime 
now occurs close to Vl = at td — 0.2. The noise cross correlation acquires an appreciable amplitude when Vl 
is increased beyond the characteristic energies associated with the double peak structure in the density of states 
(typically Vl = 0.3 and Vl = 0.7 for td = 0.2). At the location of these peaks one observes some features in the 
crossed correlation signal. For larger tunneling between the dots this monotonous decrease of the cross correlation, 
and its (larger) amplitude is even more pronounced. It can even occur for Vl < because the density of states of 
both dots acquires a peak at negative energies, thus favoring EC process as well as direct tunneling. 



V. CONCLUSION 



To summarize, wc have dcvcloppcd a framework for the description of the transport properties of a device consisting 
of a superconducting finger connected to two normal metal leads via two quantum dots adjacent to the superconductor. 
Two recent experiment o 35 i 36 measured the branching currents in such a device, and showed that non local effects were 
at play, which is consistent with Cooper pair splitting. Yet it is now well accepted that the evidence of Cooper pair 
splitting in such a device would be more robust if current-current crossed correlations were measured, which is the 
main justification for the present work. 

For pedagogical reasons, we started our analysis with the description of a device where the superconductor is 
connected to a single dot (which is in turn connected to normal metal leads), in order to derive the self- energy which 
arises from the coupling of the dot to all leads. Next, we focused on the experimental geometry of Refs. [H and l36l with 
two dots, allowing for direct tunneling between the two. The current and noise crossed correlation were computed 
using the Keldysh formalism. The branching currents were expressed in terms of the single particle Green's function 
of the dot, which is dressed by the coupling to the leads, in the same spirit as the Fisher-Lee formula^. Similarly, 
the noise (and in particular the noise crossed correlations) was expressed in terms of a two particle (dressed) Green's 
function for the dot variables. This second result corresponds to an extension of the Fisher-Lee/Landauer-Biittiker 
formula for the noise. Provided a given choice of interactions within the dot (Coulomb, electron-phonon,...) these 
formulae for the currents and current-current crossed correlation are exact if the single and two particle Green's 
functions can be computed. For our purposes, we chose to illustrate the operation of this device by ignoring such 
interactions. In this case the two particle Green's function can be decoupled as a sum of products of single particle 
Green's functions. 

The current and the noise crossed correlation were obtained numerically by solving the Dyson's equation. Given the 
complexity of the device, the number of parameters (voltages, energy levels, resonance widths) constrained us to focus 
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on few specific cases. We focused on transport for voltages and energy levels all contained within the superconducting 
gap, because this corresponds to the regime where Cooper pair splitting is understood to occur. We looked at two 
configurations for the dot energy levels: the antisymmetric configuration (the two levels have opposite energies with 
respect to the superconducting chemical potential) and the symmetric configuration (same energies for the two levels). 
A crucial concept for the understanding of the generated data for current and noise is the density of states of the dots: 
due to the proximity effect induced by the presence of the superconductor, the density of states of a given dot acquires 
some weight ("small" resonance) at the energy which is opposite to the resonant energy of the bare dot ("large" 
resonance). Electron transport is therefore favored if the electrons can choose to be transfered via either a "large" 
or a "small" resonance of these dots. We fixed the voltage of a given lead within the gap, below all resonances, and 
we varied the voltage of the other normal metal lead within the whole gap. We started the discussion assuming no 
direct coupling between the two dots. On one hand we observed that CAR process is optimized in the antisymmetric 
case, when both normal lead voltages are below the dot density of states resonances. On the other hand we observed 
that EC process is the most important one in the symmetric configuration, when the voltage which is varied is set 
above all resonances. A systematic study of the current/voltage characteristic and of the noise/voltage characteristic 
allowed to identify precisely which processes (CAR, DAR, EC) are into play when the voltage is increased. We 
argued that both the monitoring of the branching currents as well as the noise crossed correlations are necessary to 
identify these processes: for instance, DAR processes contribute to the current, but not to the noise correlation signal. 
Finally, we switched on the direct coupling between the dots and we repeated the analysis. The density of states of 
the dots undergoes strong modifications in the presence of such a coupling. This coupling has a tendency to spoil 
the positive correlations: for sufficiently strong coupling, positive correlations disappears altogether. We therefore 
provided a rather complete picture of the operation of this Cooper pair beam splitter which could be useful in future 
noise crossed correlation experiments. 

Several extensions of this work can be envisioned. First, we clearly neglected the separation between the two injection 
locations of the superconductor/dot interfaces. In all theoretical and experimental investigations, this separation r 
has to be smaller than the superconducting coherence length (the "size" of the Cooper pair in the superconductor), 
because this leads to an exponential decay of the CAR (and EC) process. In addition to this reduction, a power law 
decay with kpr (which depends on the dimensionality) is typically found theoreticall y 6 ' 20 ' 21 ' 24 and it can lead to a 



strong suppression of CAR and EC processes. Ref. [35| did not find conclusive evidence of this power law suppression. 
It argued that the segment of the nanowire which is buried below the superconductor, which acquires a minigap, 
plays the dominant role for Cooper pair splitting and should be immune to the power law decay because electrons are 
directly injected to the dots on both sides. Nevertheless, further investigations of these power law suppression with 
our setup could provide additional insight. 

Second, we have neglected the Coulomb interactions within the dot, which is justified if the resonant linewidth 
is sufficiently large (but not larger than the superconducting gap) so as to constitute open quantum dots. The 
inclusion of interactions within the dots constitutes a definite challenge, and such interactions can in practice only be 
included via approximate treatments such as mean field theory, perturbative diagrammatic resummation, or Kondo 
phenomenology. In practice, electron interactions allow to transfer electrons one by one through each dot. In our 
device, we can only argue that successive Cooper pairs emitted by the superconductor do not overlap significantly 
if the boundary between the superconductor and the dots is sufficiently opaque. In the two experiments which are 
relevant to this work, it was possible to characterize the dots by switching off the superconductivity by applying a 
magnetic field. Such analysis clearly showed the presence of Coulomb diamonds, and the boundaries between them 
allowed to identify parameters for resonant electron transfer. Here we can only argue that regardless of the physical 
origin of such resonances, they will have a characteristic location and width associated with interactions, and they can 
be effectively tuned by adjusting the gate voltages on the two dots. The scenario for adjusting them antisymmetrically 
or symmetrically to favor CAR/EC processes should then be robust. But given the fact that our intermediate results 
for the current and noise can be expressed in terms of exact one and two particle Green's functions, an extension of 
the present work including interactions using an approximate scheme is foreseeable in the future. 

Finally, we have focused solely on the zero frequency noise crossed correlation signal. Noise crossed correlations 
at finite frequencies contain additional information, such as the relevant time scale for Cooper pair splitting in this 
particular device. Several entanglement scenarios have addressed the importance of short times/high frequencies, for 
instance when performing an entanglement diagnosisj 45 ' 46 Further investigations along these lines would prove useful. 

Note added in proof. Upon completion of the manuscript, we failed to include some relevant references on master 
equations type approaches which are applied to study transport in the entangler of Ref. 0. Refs. H3 and 4§ derived 
the quantum master equation for the entangler and computed the branching currents in this device. Ref. 48| focused 
on the experimentally measurable quantities in the light of the recent experiments Refs. l35l and[36l Ref. |4S computed 
the zero frequency noise crossed correlator in order to achieve a Bell inequality test within this approach. 
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